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The standard calculations of the ground-state energy of a homogeneous Bose gas rely on 
approximations which are physically reasonable but difficult to control. Lieb and Yngvason 
[Phys. Rev. Lett. 80, 2504 (1998)] have proved rigorously that the commonly accepted leading 
order term of the ground state energy is correct in the zero-density- limit. Here, strong indications 
are given that also the next to leading term is correct. It is shown that the first terms obtained in 
a perturbative treatment provide contributions which are lost in the Bogoliubov approach. 
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As dilute weakly interacting Bose gases are experimentally realisable since 1995 B , there is a renewed interest 
in principal results for the ground state energy derived some 50 years ago In recent years, several 

subtle issues have been clarified [ill IT^ IT^ II^ . Consider a weakly interacting translationally invariant Bose gas in 
the thermodynamic limit with non-negative two-particle interaction potentials excluding bound states of two or more 
particles. In momentum representation, the A^-particle Hamiltonian in second quantisation reads (15» | 

H = ^^aia^ + \y {kk2\U\k^k)a\a\a^a^ , (1) 
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where the matrix elements 



{hk2\U\hh) = 5-^^-^_-^^-^[/(fc2-fc4) (2) 
are given by the Fourier transform of the potential 

U{k) = ^ j C/(r)e'*^-^-dV . (3) 



V 
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This approach is valid both for soft potentials (non-pseudo-potentials for which the Fourier transform exists) and the 
pseudo-potential usually used to model hard sphere interaction 

, , Airah^ ^, , d 

Uoif) = Sir)—r , 4 

TO or 

where a is the hard sphere diameter. 

For the pseudo-potential I^J, neglecting terms of the Hamiltonian which are believed to be small leads to the 
Lee- Yang formula P, (Bl for the ground-state energy per particle eo = Eq /N 



2Trah'^ 
eo = n 



15V^ 



(5) 



where n = N/V is the density. The o na^^ includes a term proportional to na"^ In (na'^) as well as higher order 

■ terms (see ref. and references therein) which are very small for realistic experimental values of na^ . According 
. ^ . to the so-called "Landau postulate" , eq. © is the ground state energy for any non-negative interaction potential if 
the hard sphere diameter a is replaced by the s-wave scattering length [T^IT^. This is supported by the fact that 
JH . known potential dependent corrections to eo are given by Cna^2TTah^n/m [8|. The constant C has only recently been 
■ - - ■ calculated explicitly by Braaten et al. in terms of a quantity defined by the 3^3 scattering amplitude 0, 0| . 

Lieb et al. [lUlSj have calculated rigorous bounds to prove that in the limit of vanishing density n, the leading order 
term of the ground state of the complete Hamiltonian is indeed given by the leading term of the Lee- Yang formula: 

1 - C(na3)Vi7 < < 1 + ii.3(na3)i/3 ,00. (6) 
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However, despite many efforts, the next to leading term could not be proved rigorously; neither its sign nor the 
exponent of na^ has been established beyond doubt. 

A standard textbook derivation of equation lO uses the Bogoliubov approximation to replace the quartic Hamilto- 
nian J^l by a quadratic Hamiltonian which can be diagonalised exactly with the help of the Bogoliubov transformation. 
Repeating the same analysis for soft potentials leads to the formula of Brucckncr and Sawada j6i, illl IT^ 

(0) 2T:h'^{ao + ai) 2nh'^ao 128 / T 

where aq and ai are the first and second Born approximation of the scattering length. In general, the s-wave scattering 
length a is not identical with the sum of the first two Born approximations ao + ai. Thus, the terms of the Hamiltonian 
neglected in the Bogoliubov approximation have to contribute in leading order in the density to the ground state energy 
per particle. 

In this letter we use quantum mechanical perturbation theory to show how the missing terms of the Born series 
can be recovered. Furthermore, for very soft potentials, already first order perturbation theory indicates that the 
next to leading term of the ground state energy of the full Hamiltonian is also given by the term of the Lee- Yang 
formula. If one accepts the Landau postulate, any experimentally realistic interaction potential can be replaced by a 
soft potential with a rapidly converging Born series. However, our approach is not restricted to such potentials; an 
extension to non-regular potentials like the pseudo-potential Q will be given in a subsequent paper. 

We also improve the upper bound © for the ground state energy. Lieb et al. used a product-ansatz for the 
iV-particle wave function to calculate their upper bound according to the variational theorem. Here, we use the 
Bogoliubov ground state (which is the exact quantum-mechanical ground state of an approximate Hamiltonian) to 
calculate an upper bound on the ground state energy of the complete Hamiltonian. 

Although we do not use the Born approximation for the scattering length at any point in our calculations, several 
expressions can be identified with terms of the Born series which are derived in a first step. Before we can calculate 
the upper bound, we describe the properties of the Bogoliubov ground state for soft potentials. 

I. BORN SERIES 

Often, the scattering length is approximated by the first term of the Born series 

The next two terms of the Born series are given by |^ 

\r- fi 



, 3 

a2 



\r -ri\ \J \ri - r2\ 
higher order terms can be calculated analogously. The identity 



dV ; (10) 



1 4-Tr pik-{r-ri) 

Y.—^+^^^") (11) 
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is useful to rewrite the terms of the Born series to expressions in momentum space. Using eq. 0, we obtain both the 
well known formula (see e.g. ref. |17|') 



MO 



and 

V /TO\3 v-^ U{ki)U{k2)U{ki~k2 



(my ^ U fci U fc2 fei-fc2 
-^-tAy^)X. ^2 ■ (13) 
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II. 



BOGOLIUBOV GROUND STATE 



For any (normalised) A^-particle wave-function |'0o)i the variational theorem yields an upper bound on the ground 
state energy per particle eg of the complete Hamiltonian ^ 



eo < jjiiJoim^Po) 



(14) 



Here, we use the Bogoliubov ground state as our \tlJo). The idea of the Bogoliubov transformion is to express the 



operators al and a- which create or annihilate particles with momentum k by quasi-particle operators bl and 6-: 



with 



L fcl ' k2' l^lM ' ^ fcl ' fe2 ^ 



(15) 
(16) 

(17) 



Without the unnecessary approximation U{k) ~ U{d) (wh ich leads to unphysical divergences in both the Born series 
and in the calculations for the ground state energy |l7t Il9||). the (real) transformation parameters for isotropic 
interaction potentials (for which U{k) = U{—k)) are defined via: 



s{k) 



= 1 



k 

2m 



L-k 



1-Ll 

k 



e{k) 



NU{k) 
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h'k 
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NU{k)^ 



1/2 



(18) 
(19) 
(20) 



In the Bogoliubov ground state 1-00) on average No particles are in the ground state characterised by k 
macroscopically occupied ground state allows a grand-canonical description of the excited states. One has 



(^0|St6g|^0) 



= No, 
= No, 

= for k^O , 



0. The 

(21) 
(22) 
(23) 



where the first equation defines A'o whereas the second is only true in the thermodynamic limit when No 

is neghgible). 



as — > oo (such that the error in the approximation y^No + 1 ~ VNo is negligible). It relies on the fact that 
the number of particles in the ground state fluctuates for a weakly interacting Bose gas even at zero temperatures. 
The contribution to the total wave function of wave functions with exactly No + v or exactly A'o -1-;^ -1-1 particles are 
practically the same. However, we would like to stress at this point that the calculation of an upper bound does not 
use special assumptions about the ground state — it is of course possible to construct a wave function which satisfies 
eqs. to (1^ . A wave function with a phase factor in eq. or even with (?/jo|6-6-|V'o) = would lead to a higher 



upper bound and thus can be discarded here. 



III. CALCULATING THE UPPER BOUND 

Let us now write the full Hamiltonian ^ in the form 



H — Hi + Hii + Hill + Hiy, 



(24) 
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where the first part includes both the kinetic energy and those terms with all four ki = 0, the second those contributions 
with two momenta equal to zero. In the third part, Hm, only one ki vanishes whereas in the last all are different 
from zero: 



1 fe2L2 

Hi = -C/(d)atata-a- + V ata- (25) 

fc 

I E {umial^a^a^ + U{k)ala^^.^a_.^ + C/(fc)atat a^^a^ 

fc#0 

+ U(k)aiaia^a^ + U(0)aiala^a^ + U(0)alaia^a^] 

^ ' k k ^ ' k k ^ ' k k) 

1 y {kik2\U\hki)ai 4 a- a- . (26) 

2 ^-^ ' ' ' ' ki k2 k3 ki 



2 

{fc,#0} 



We have 



{ijo\HuM=0 (27) 



as the total number of creation and annihilation operators for k ^ is odd. Thus, the precise equation for Hm is not 
relevant here. 

All sums are expanded in powers of n with the following procedure: if (^) • ■ • finite and the next derivative 
with respect to n diverges, the expression J2 ■ ■ • is evaluated. In the thermodynamic limit, the resulting 

integrals can be evaluated for low densities after the substitution fi^k^ = 2mNU{6)y'^ 0, 0|. To calculate the 
depletion N — Nq, we also apply the fact that a complete set of eigenfunctions was used to obtain the second 
quantised Hamiltonian in the form Thus, the conservation of the total number of particles N can be expressed 

as 

N = (V'oiata^^l^o) + E(V'ol4«fel'^o) • (28) 
Because of a^j = 6g and {ipo\bVb -.lipo) = Nq , this leads to 



^ "^0 

iV- 



No = yvl^N^^,. (29) 



k^O 



The fact that in this formula for the depiction, oq rather than a appears, is due to the approximate ground state 
wave- function used at this point. 

Even for the anomalous condensate fluctuations of a homogeneous Bose gas [22, 113, 0| we can re- 
place {^po\alaia^a^\'ip())/N by {ipolala^lipo)'^ / N for large N. We get including the order O (^ny^nalj 

(V'ol-ffilV'o) _ Stt^i^oo^ 2TTh'^ao _ 16 64 



- -n n n—=Jna-^ n—^Ja-^n (30) 

and 

(V-oliJiil^o) „27rft2ai f2Trh^ao 16 2Trn^ai 16 , 2Trh^ao 32 \ 

= 2 n -I- n—^ H n—^ \ af.n . (31) 

N m \ m y/TT m 3y/n m 3y/n y v " 

The only non-zero contributions to jj{ipo\HiY\ipo) come from expressions where at least two ki have the same modulus. 
We find: 

(V'oliJivlV'o) _ 1 V- 1 p.,^ _ 2TTh^ 16 2nfi- 



1 V- 1 Tr.r 27rri2 16 2Trh'^ I — - 



N N ^ 2 

k,k'^d 

as the sums ;^ X^fc /p^^o I^i^P^^*^) ^^'^ w'^k k'^o^^'^^'^'^^ ''^ ^^'^ °^ order n^ajj. Thus, including orders 
of O (^ny^nalj S^^ ^he upper bound {^^:Q\H\'lpo) /N: 



2T:h^(ao + ai + a2) 2T:h^ao 128 / — ~ 2-Kh^ai 2,2 , ^ , , 

Co < nH n——=\ na'A^ n--—=\ na9,, (33) 

m m 15v7r V " rn 3^/71 ^ 
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which hes above the Lee- Yang formula (for non-negative potentials, the Born series is an alternating series with a < 
Qq + ai + ■ Had we chosen a wave- function with {ipQlbJj^lipo) = , the first three terms in eq. H31|l would be missing. 
The leading order in eq. would then read 27r/i^(ao — ai + a2)n/m which because of oi < lies above our bound. 
For physically reasonable densities of up to 128(na'^)^/^/(15-y/7r) « 0.03 |0 and a potential such that the Born series 
converges fast enough, we have (ao + ai + a2 — a) / a <^ 11.3(na'^)^/'' « 0.4. Thus, with respect to the exponent of the 
next to leading term, our bound constitutes a significant improvement of the upper bound © . 



IV. GROUND STATE ENERGY BEYOND THE BOGOLIUBOV APPROXIMATION 



We note that the Bogoliubov transformation exactly diagonalises the Hamiltonian Hq ~ Hi + Hu (see e.g. ref. 

AT ^ fc2 1,2 ivr , ^ ^ ~i 

k k^O 

fe/0 

with quasi-particle energies e{k) given by eq. H2()|l and a ground state energy per particle eg'^^ given by the Brueckner- 
Sawada- formula Q). We now employ Hq as the starting point for perturbation theory, 

H = Ho+[h- Ho) (35) 

and again assume an inter-particle interaction potential such that the Born series rapidly converges. Then already 
the first order correction to the ground state energy (c/. eq. (|33|l ') 

(1) _ 1 , , , ^\ ^\ : , ^ '2Trh'^a2 2TTh^ai 32 



A* mm 3V7r 

is small. 

Analogously, it can be shown that in second order perturbation theory we have 

(2) _ 2TTh?as, 



e);> = n + O [n^nalj (36) 

as the leading order term. Thus, the terms of the Born series which are missing in cq. ^ are recovered, for low 
densities, by Rayleigh-Schrodinger perturbation theory, i.e. they stem from contributions to the Hamiltonian which 
are neglected in the usual Bogoliubov approximation. Analogously, the ground state wave function can be modified 
to obtain an improved upper boimd. 



V. CONCLUSION 



Using a quantum-mechanical variational ansatz, we have derived an upper bound H33|l on the ground state energy 
which lies below the bound ^ for physically reasonable densities if the Born series converges fast enough. We have 
shown how our approach could in principle be used to further improve this bound. However, more important than 
an ever improved upper bound is the fact that that there are strong indications that a perturbative treatment of the 
ground state energy starting with the decomposition (|35|l . actually converges towards the Lee- Yang formula ^ (if 

both the Born series and a'^/^ = ao^^(l + ai/ao + 02/00 -I- . . .)^^^ — ao^^(l -I- |ai/ao -f- . . .) converge). 
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